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The SU{2)l x SU{2)ii chiral quark model consisting of the (cr, tt) meson multiplet and the con- 
stituent quarks propagating on the homogeneous background of a temporal gauge field is solved at 
finite temperature and quark baryon chemical potential fiq using an expansion in the number of 
flavors Nf , both in the chiral limit and for the physical value of the pion mass. Keeping the fermion 
propagator at its tree-level, several approximations to the pion propagator are investigated. These 
approximations correspond to different partial resummations of the perturbative series. Comparing 
their solution with a diagrammatically formulated resummation relying on a strict large- A^/ expan- 
sion of the perturbative series one concludes that only when the local part of the approximated pion 
propagator resums infinitely many orders in 1/Nf of fermionic contributions a sufficiently rapid 
crossover transition a.t fiq — is achieved allowing for the existence of a tricritical point or a critical 
end point in the fiq —T phase diagram. The renormalization and the possibility of determining the 
counterterms in the resummation provided by a strict large- A^/ expansion are investigated. 
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I. INTRODUCTION 
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The low-energy effective models of the QCD, such as the Nambu-Jona-Lasinio (NJL) model [1] and the linear sigma 
■ model (also called chiral quark or quark-meson (QM) model) [2] are based on the global chiral symmetry of the QCD. 
. ] They proved to be very useful in qualitative understanding of many aspects related to the spontaneous breaking of 
, the chiral symmetry and its restoration at finite temperature and density, but share as a major drawback the lack of 
I ' the confinement property. As a consequence of the absence of gluonic effective degrees of freedom and due to the lack 
, of color clustering [3] there are unsuppressed contributions of constituent quarks in the low-temperature phase. Both 
' features, which are in fact related, alter the reliability of the quantitative thermodynamic predictions of these models, 
' such as the equation of state or the location of the critical end point (CEP) in the fiq — T phase diagram. 

Since the QCD phase transition involves both the restoration of chiral symmetry measured by the evaporation of the 
chiral condensate (i'lp) and the liberation of quarks and gluons encoded in the change of the Polyakov loop $, much 
T-H . effort was devoted to understand the relation between the chiral and deconfinement phase transitions. An argument 
' by Casher [4] states that in the vacuum, that is at /i^ = T = 0, confinement implies the breaking of the chiral 
.symmetry. The connection between {'iptp') and $ is revealed by the spectral density of the Dirac operator p(A). In 
the vacuum the Banks-Casher relation p(0) = {'tpip) /i" [5] states that the spectral density of the Dirac operator in 
. the deep infrared is proportional to the quark condensate. Finite temperature lattice studies at fiq — show that the 
' connection between chiral symmetry and confinement suggested by Casher's argument holds. The infrared part of p(A) 
undergoes a pronounced change as one crosses from the confined to the deconfined phase [6, 7]. As shown recently in 
[8], the phase of the Polyakov loop $, which can be expressed as a spectral sum of eigenvalues and eigenvectors of the 
Dirac operator with different boundary conditions, receives its main contribution from the infrared end of p(A). It is 
generally true for both Nc — 2 and Nc — 3 that at high temperature the fermion determinant favors the sector where 
the Polyakov loop lies along the positive real axis [6, 7]. For this type of configuration in which the phase of $ vanishes 
the chiral symmetry is restored, because a sizable gap develops in the spectral density of the Dirac operator which 
implies in view of the Banks-Casher relation {■ip4'} = 0- For configurations in which the phase of the Polyakov loop is 
not vanishing, the chiral symmetry is not restored, as observed in the lattice study of the quenched QCD [9] and also 
by using a random matrix model calculation [10]. 

Casher's argument suggests that the temperature Td for the deconfinement phase transition is somewhat lower than 
the restoration temperature of the chiral symmetry. As explained in [11] at finite density Casher's argument could 
fail, so that it does not contradict the existence of a dense phase in which at a given temperature chiral symmetry 
is restored while quarks remains confined. Such a phase can exist inside the so-called quarkyonic phase, which was 
suggested as a new phase of the QCD at finite temperature and density, based on its existence within a large-TVc 
analysis [12]. 
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In the Polyakov-loop extended NJL model (PNJL) where the couphng of the Polyakov loop to the quark sector is 
achieved by the propagation of the quarks on a constant temporal gauge field background, the simultaneous crossover- 
type transition of deconfinement and chiral restoration was obtained [13]. As shown in [14, 15] this model is able 
to reproduce the main features of the quenched lattice result of [9]. The phase transition was recently intensively 
investigated in the PNJL model with two [16-20] and three flavors [21-23], also in the nonlocal formulation of the 
model [24, 25] . The interplay between chiral and deconfinement transitions was investigated in the PNJL model using 
large- iVc expansion in [26]. 

By coupling the Polyakov loop to the quark degrees of freedom of the QM model the thermodynamics of the resulting 
Polyakov quark- meson model (PQM) was studied for two [27-32] and three quark flavors [33-35]. The effect of a strong 
magnetic field expected to be generated in the LHC in noncentral high-energy heavy ion collisions on the chiral and 
deconfining phase transitions was studied recently in [36] within the PQM model with two flavors. The possibility of 
coupling the Polyakov loop to meson models without quarks was considered in [37, 38]. 

The coupling of the Polyakov loop to the chiral effective models mimics the effect of confinement by statistically 
suppressing at low temperature the contribution of one- and two-quark states relative to the three-quark states. This 
feature makes the Polyakov-loop extended effective models more appropriate for the description of the low-temperature 
phase and for quantitative comparison with the thermodynamic observables on the lattice [27, 33, 39]. Better agreement 
is expected up to T ~ (1.5 — 2)Tc above which the transverse gluonic degrees of freedom dominate in thermodynamic 
quantities, such as the pressure, over the longitudinal ones represented by the Polyakov loop. 

Despite this success, one should keep in mind that the solution of the Polyakov-loop extended effective models is 
mainly obtained in the lowest one-loop (trace-log) order of the fermionic sector, hence studying their stability against 
inclusion of higher loops would be certainly of interest. Some approximations to the PQM model [27, 28, 33 35] neglect 
the fermionic vacuum fluctuations and by treating the mesonic potential at tree level, completely disregard quantum 
effects in the mesonic sector. The effect of including the quantum fluctuation in the PQM model was recently studied 
in [29, 30, 32] using functional renormalization group methods. 

In this work we would like to address two questions using large- A^y approximation to the SU{2)l x SU{2)fi ~ 0(4) 
model. The first one is to what extent the inclusion of the Polyakov loop modifies the Hq — T phase diagram obtained 
previously in [40] in the chiral limit of the two flavor QM using the large- iV/ approximation. The second one concerns 
the effect of different partial resummations on the quantitative results. To this end several approximate resummations 
of the perturbative series will be investigated and the obtained results compared. 

The paper is organized as follows. In Sec. II we overview some basic facts about the Polyakov loop, including different 
forms of the effective potential and we introduce and parametrize the PQM model, presenting also the approximations 
exploited for its solution. The renormalization of the model and the determination of the countertcrins is discussed 
in Sec. III. In Sec. IV we present the numerical results on the Hq — T phase diagram obtained in the chiral limit and 
for the physical value of the pion mass by using different forms of the Polyakov-loop effective potential and various 
approximations to the resummed pion propagator. Section V is devoted to discussion and summary. 



II. THE PQM MODEL WITHIN A LARGE-iV/ APPROXIMATION 



A. The Polycikov loop as an order parameter 



We shortly review a few well known facts about the Polyakov loop incorporated as a new effective degree of freedom in 
the chiral quark model. This is usually done by considering the propagation of quarks on the homogeneous background 
of a temporal gauge field Ao{x). At finite temperature T = 1//3, after analytical continuation to imaginary time t — )■ ir, 
— > 1^4, the temporal component of the Euclidean gauge field A/^ enters in the definition of the Polyakov-loop 
operator (path ordered Wilson line in temporal direction) L(x) and its Hermitian (charge) conjugate L'^{x) 



L{x) = Vexp 



dTA4{T,x) 



L^{x) = Vexp 



dTAl{T,x) 



(1) 



which are matrices in the fundamental representation of the SU{Nc) color gauge group {Nc = 3). In the so-called 
Polyakov gauge, the temporal component of the gauge field is time independent and can be gauge rotated to a diagonal 
form in the color space ^4,d(^) = <f>3{x)X3 + ^s{x)X8 [41-43], where A3, Ag are the two diagonal Gell-Mann matrices. 
Then the Polyakov-loop operator simplifies 



L{x) = diag(e' 



i/30+(x) _i/30-(x) .-i/3(0+(x)-|-0_(x)) 



(2) 



where 4'±{x) = ±(f>3{x) + (f>six)/y/3, with a similar form for the conjugate L^{x). 

Topologically nontrivial gauge transformations U{t, x) € SU{Nc) that are periodic up to a twist, that is U{t+(3, x) = 
zU{t,x), were introduced in [44]. where z is an element of the center of the SU{Nc) group which is isomorphic with 
Zat^ = {z\z = cxp {27mi / Nc) , n = 0, 1, . . . , A^c ^ 1}, the cyclic group of order Nc- Under such transformations the 
color trace of the Polyakov-loop operator and its conjugate, that is l{x) = tTcL{x)/Nc and /^(x) = trcL^(a;)/A^c, are 
transformed by an element of the center: l{x) zl{x),l^x) — >■ z*l^x). Consequently, in the pure gauge theory. 
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which has an exact Z;Vc global symmetry, the thermal expectation value of the traced Polyakov-loop operator and its 
conjugate, 

m = i- (treL(f)) , ^x) = (tr,it(f)) , (3) 

are order parameters for the center symmetry and must vanish if the symmetry is unbroken. However, the Polyakov loop 
(t>(.T) and its conjugate can acquire a nonvanishing value, signaling the spontaneous breaking of tlie Zjv^ symmetry. 
These complex quantities can be regarded as order parameters of the deconfinement phase transition, because the free 
energy of a heavy (static) quark-antiquark pair with spatial separation f = x — y is related to the expectation value 
of the correlator of the traced Polyakov-loop operator for which cluster decomposition is expected to hold at infinite 
separation 

exp [-/3F,,(r)] = {L{x)L\x)) ^ (4) 

When i>(y) = then Fqq{r) — >■ cxd, and when $(y) ^ then Fqq(r) is finite, which means confinement and 

deconfinement, respectively [45, 46]. 

In the presence of dynamical fermions the Z^v^ symmetry is not exact anymore. Nevertheless, the Polyakov loop gives 
through its distribution information about the confinement (low T) or deconfinement phase (high T) of the system both 
in a canonical or grand-canonical formulation of the QCD [47, 48]. Since its absolute value can be related to the free 
energy difference between two systems, one containing the quark-antiquark source pair and the other not containing 
it, by renormalizing the free energy a renormalized Polyakov loop can be defined [49], which provides information on 
the temperature of the deconfinement phase transition. 



B. The mean-field Polyakov-loop potentials 

In the mean-field approximation <I>(x) and $(x) are replaced by .^-independent constant fields which satisfy |$| = |<l| 
at vanishing chemical potential. We review here several forms and some basic features of the mean-field effective 
potential for the Polyakov loop frequently used in the literature. This effective potential will be incorporated in the 
thermodynamic potential of the PQM model. The simplest effective potential is of a Landau type, constructed with 
terms consistent with the Z3 symmetry [50]: 

^^ZYpo,y($, = - |($3 + $3) + |(^5)2 ^ (5) 

where the temperature-dependent coefficient which makes spontaneous symmetry breaking possible is 



62(T)=ao + ai(^^j+a2(^^j • (6) 

To is the transition temperature of the confinement/deconfinement phase transition, in the pure gauge theory Tq = 
270 MeV. The parameters ai,i = 0, ... ,3 and 63,64 determined in [51] reproduce the data measured in pure SU{3) 

lattice gauge theory [52] for pressure, and entropy and energy densities. The minimum of this potential is at $ = 
for low temperature and $ — > 1 for T — > 00 in accordance with the definitions (1) and (3). However, when using this 
potential in either the PNJL or the PQM models the minimum of the resulting thermodynamic potential is at $ > 1 
for T ^ 00 and also leads to negative susceptibilities [24] . 

An effective potential for the Polyakov loop inspired by a strong-coupling expansion of the lattice QCD action was 
derived in [13, 53]. Using the part coming from the SU{3) Haar measure of group integration an effective potential was 
constructed in [54, 55] 

/3'*Wiog(*,3>) = -ia(T)$i>-h6(r)ln l-6$l'-F4($3 + $3^ _3^^^^2 ^ 
with the temperature-dependent coefficients 



a(T) = ao + ai(^^j+a2(^^j , 6(r) = 63 (^^ j . (8) 

The parameters ai^i = 0,1,2 and 63 determined in [54] reproduce the thermodynamic quantities in the pure SU{i) 
gauge theory measured on the lattice. Tlie use of this effective potential cures the problem with negative susceptibilities 
[24] and since the logarithm in U\og{^, $) diverges as $ ^ 1 it will also guarantee that $ ^ 1 for T ^ 00. 
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A third effective potential is the one determined in Refs. [13, 53]: 

/3WFuku($, = -b [54e-"/^$l> + hi [l - 6$$ + 4($'' + ¥') - 3($l>)2]| , (9) 

where a controls the temperature of the dcconfincmcnt phase transition in pure gauge theory, while b controls the 
weight of gluonic effects in the transition. In this case, the parameters a = 664 MeV and b — (196.2MeV)'' are obtained 
from the requirement of having a first order transition at about T = 270 MeV [27, 56] . 

It was shown in [56] that there is little difference in the pressure calculated from the three effective potentials for 
the Polyakov loop in their validity region up to T ~ (1.5 — 2)rc- The presence of dynamical quarks influences an 
effective treatment based on the Polyakov loop which in this case is not an exact order parameter. Defining effective 
Polyakov- loop potentials for nonvanishing chemical potential when |<I>| ^ |l>| is somewhat ambiguous [27]. In the 
present analysis we will use at /ig ^ the Z3 symmetric Polyakov-loop potentials given above. The effect of the 
dynamical quarks was modeled by considering the Nj and fig dependence of the To parameter of the Polyakov-loop 
effective potential. Using renormalization group arguments this dependence was parametrized in [27] to be of the form 
To{iiq,Nj) = Trexp(-l/(ao6(/x,))). The parameters were chosen to have To(/iq = 0,iV/ = 2) = 208.64 MeV. When 
using the polynomial and logarithmic effective potential for the Polyakov loop given in (5) and (7) we will consider 
in addition to Tq = 270 MeV two more cases, one with a constant value Tq = 208 MeV and the other with the 
above-mentioned /i^-dependent Tq taken at Nf = 2. 



C. Constructing the grand potential of the PQM in a "^-derivable" approximation 

The Lagrangian of the SU{2)l x SU{2)ii chiral quark model [2] is written in the usual matrix form [57, 58] by 
introducing two Nf x Nf matrices 



M{x) 



ll + iT"7r"(a;), M^^x) 



l + ^75TV(x), 



(10) 



in terms of which one has 



C = tTf [d^M'^d'^M - rn^M^M] - — [try {MH'l)f + Nfha + ip {ij^'Df, - gM^) ij) + c.t. 



(11) 



where in the mesonic part we have introduced an explicit symmetry breaking term and "c.t." stands for counterterms. 
After performing the trace, one can see that without the fermionic term the Lagrangian (11) is that of the 0{N) model, 
which describes the system of sigma and A'^ — 1 pion fields and is appropriately parametrized for a large- A^ treatment 
[59]. Vanishing background is considered for the spatial components of the gauge field and a constant mean- field Aq 
for the temporal component, so that the covariant derivative is = — IS^qAq. The trace in (11) is to be taken 
in the flavor space and to simplify notations the flavor, color, and Dirac indices of the fermionic fields xj), are not 
indicated. The SU{Nf) generators in the fundamental representation (a = 1, . . . , N'j — 1) satisfy the normalization 

condition tr/ (T^T**) = S"'^ /2. Some rescaling with Nf = y/N was done and since in the mesonic sector we only want to 
increase the number of pions we do not introduce another invariant, tr/ [(M^'M) ] , which for Nf > 2 is independent 
of [tr/(M'l'M)]^. For a recent treatment of the U{Nf)L x U{Nf)R meson model having both invariants see [60]. 

The constituent quarks become massive only after spontaneous/explicit symmetry breaking. In (11) the sigma field 
is shifted by the vacuum expectation value v as a ^ v^/N + a, where on the right-hand side of the arrow a denotes 
the fluctuating part of the original field. Then evaluating the trace, one obtains 



jC = -N 



'^4,122 u 

— V H — m V — nv 
24 2 

1 



- /a 



V +m V — h 



+Vi[(i9'' + ^''°Ao)7^-mJV- 



-m 

9_ 

A 



-.up 



^sfN"" 24N 
Tp(a + i^/Wf^^T^-TT 



'ij; + c.t.. 



(12) 



where = a"^ + -jf^ . Here a rcscaled Yukawa coupling g = g^/N/{2Nf) was introduced in order to assure the finiteness 
of the tree-level quark mass niq = gv in the N ^ 00 limit. In this limit, due to the Nf scaling of the vacuum expectation 
value of the sigma field in (10) and of the coupling A in (11) the tree-level sigma and pion masses m^g = m'^ + At;^/2 
and m^Q = + are also finite. 

After continuation to imaginary time, the grand partition function Z and the grand potential i^{T, hb) oi the spatially 
uniform system defined by (12) are introduced as follows: 



Z = tr{ exp [ - PiHoiAi) + Hint - Mb<3b)] } = e" 



(13) 



where hb is the baryon chemical potentiaL 7Ji„t is the interacting part of the Hamiltonian constructed from (12). 
With the help of Hq, the quadratic part of the Hamiltonian at vanishing A4, one can define the A4-dependent free 
Hamiltonian iJo(^4), which for two quark flavors u and d reads as 

H^iA.) =Ho + J + , (14) 

where A4 = diag(0+ , 0_ , — ((/>+ + (p- j) is diagonal in color space and i,j denotes color indices. In (13) Qb is the 

3 

conserved baryon charge which can be written in terms of the particle number operators as Qb = ^ J2 ^g.ii ^ith 

i=l 

Nq^i = Nu^i + Nd.i — Nu^i — and e.g. Nu,i = j d?x{u\ui + d\di). Then combining Hq{A4) and ^bQb one 
can see that the effect of fermions propagating on the constant background A^, diagonal in the color space, is like 
having imaginary chemical potential for color. Following Ref. [61] one introduces color-dependent chemical potential 
for fermions 

Hi,2 = l^q - i(p±, Hz = Hq + i{(j)+ + (!>-), (15) 

3 

where = is the quark baryon chemical potential. Then, introducing the notation % = Hq—^ HiNq^i, one can 

i=l 

write Z as a path integral over the fields, generically denoted by ^ 



J [25*]|e-^«Pexp[-^ (irirint(r)]| 



(16) 



where flo is the grand potential of the unperturbed system with fermions having color-dependent chemical potential 

^0 =y [v^]e-I^H. (17) 



In the onc-particle irreducible (IPI) formalism the grand potential is given by the sum of the grand potential of 
the unperturbed system and of perturbative quantum corrections represented by closed loops constructed with the 
tree-level (free) propagators. In the "^-derivable" approximation of Ref. [62], also called two-particle irreducible (2PI) 
approximation, the grand potential is a functional of the full propagators and field expectation values, and is of the 
following form: 

m[G^, G,, G, V, = f/($, $) + ^m^v^ + 7vA^4 _ ^/^^ _ _ ^ [lnG-\k) + D-\k)G^{k)] 

f [\nG-\k)+D-\k)G^ik)] +VNitvD,c I [\nG-\k) + D-\k)G{k)\ 
2 Jk Jk 

+r2Pi[G,,G,,G,«,$,$] + c.t., (18) 

where the trace is taken in Dirac and color space. C/($, $) is a particular version of the effective Polyakov-loop potential 
reviewed in Sec. II B; the tree-level propagators of the pion, sigma, and constituent quark fields are 

iD-\k) = e - mlo, iD-\k) = e - mlo, iD-\k) = ^ - m^, (19) 

while G,r, Gcr, and G are the respective full propagators. r2Pi[G^, G^-, G, $, "1] denotes the set of 2PI skeleton 
diagrams constructed with full propagators, which to 0{\/^fN) accuracy is given by 



r2Pi[G^, G,, G, V, $] = TV A ^ G^{k)^ +^j^G^{k) G^ip) - n(A;) -'-J^\n(l- ^Il{k) 



AH(fc)/6 

-s/N^liirD,c I f-f5G{k)^sG{k+p)G„{p) + -^itTD,c f f G{k)G{k + p)G,{p) , (20) 
^ JkJp 2\/N JkJp 

where the notation n(fc) = —i GTr{p)GTr{k +p) was introduced. The mesonic part of r2pi contains the 2PI diagrams 

Jp 

of the 0{N) model as given in Eq. (2.13) and Figs. 2 and 4 of [63] and also in Eq. (48) and Fig. 2 of [64]. We sec that 
the contribution of the fermions which goes with fractional powers of N intercalates between the leading order (LO) 
and next-to-leading order (NLO) contributions of the pions, which go with integer powers of N. This can be also seen 
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by comparing the expression of the pion propagator derived in (29) with Eqs. (53) and (54) of [64]. 

We use the imaginary time formahsm of the finite temperature field theory in which the four-momentum is A; = 
{iiUn, k). The Matsubara frequencies are a;„ = 2-jTnT for bosons while for fermions they depend also on the color due to 
the color-dcipcndcnt chemical potential /ij introduced in (15) and are given by u>n = (2n + l)nT — ijii. The meaning of 
the integration symbol in (20) is 



(21) 

(27r)3 ^ ' 



The dependence on $ and $ of the fcrmionic trace-log term in the grand potential VL and of the quark-pion setting- 
sun in r2Pi results from the fact that, as shown in the Appendix between Eqs. (A.l) and (A. 5), as well as between 
Eqs. (A. 35) and (A. 37), respectively, after performing the Matsubara sum the color trace can be expressed in closed 
form in terms of the mean-field (.x-indcpcndcnt) Polyakov loop $ and its conjugate $. 

What we evaluate in this work is not the grand potential (20), but rather its derivatives, that is the equations for 
the two-point functions and the field equations, which are given by the stationary conditions 

5VL _ 50. _ 50 _5Q. _5Q. _5Q. _ 

6G~6G^~5G^~1^~6^~M~ ' ^ ' 

In each of these equations we will keep the contribution of the fermions only at the leading order in the large- AT/ 

expansion. The LO contribution of the fermions is 0{^/N) in the field equations of $ and 0(1) in the equation for 
the fermion propagator G, and ©(l/v^iV) in the remaining equations, that is the field equation of v, and the equations 
of Gtt and Ga- 
it is easy to sec that the third and fourth terms on the right-hand side of (20) do not contribute to any of the 
equations at the order of interest, and that the second term contributes only in the equation for the sigma propagator 

In fact, the equation for G^ decouples, in the sense that G^ will not appear in any of the remaining five equations. 
Nevertheless, it plays an important role in the parametrization of the model, as will be shown in Sec. II E. 



D. Approximations made to solve the model 



In this work we use some approximations to solve the set of coupled equations coming from (22). 

1. As a first approximation we disregard the self-consistent equation for the fermions arising from 5SI/5G = 0, that 



IS 



iG-\k) = iD-\k) ^ tg^ / 75G(p)75G^(p-A;) + c.t., 



(24) 



and simply use in the remaining five equations the tree-level fermion propagator D{k). A study based on the solution 
of the self-consistent equation for the fermion propagator is left for a forthcoming publication. 

Within this approximation the field equation for v, hereinafter called equation of state (EoS), and the pion propagator 
simplify considerably. The contribution of the last but one term of (20) to the pion propagator breaks up upon working 
out the Dirac structure into the linear combination of a fermionic tadpole T(mq) and a bubble integral I{p;mq). 
Introducing the propagator 



these integrals are defined as 



N 

7(p;m,) = — ^ -i Do{q)Do{q+p) 



i=l •'9 

In terms of these integrals which are evaluated in the Appendix between Eqs. (A. 11) and (A. 18) one obtains: 

h] 



= Nv 



m 



c.t. , 



(25) 

(26) 
(27) 

(28) 
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iG-\k) =e -m^ -■^(v'' + ( G^{k)] + ^^f{m,) - ^^k'^Iip; m,) + c.t. . (29) 

One can see that the Goldstone theorem is fulfilled, since using the EoS in the equation for the pion propagator one 
obtains iG~^{k = 0) = —h/v. This is only accidental because the Ward identity relating the inverse fermion propagator 
and the proper vertex r^„^^ = S^T/S-tjjSipS-jT'' (see e.g. Eq. (13.102) of [65]) 



-'-Ta{j5,iG-\p)} =v^ ^T^.^^{0,p,-p), (30) 

is satisfied only with tree-level propagators and vertices. The relation above is violated at any order of the perturbation 
theory in the large- iV/ approximation, since in view of (24) the corrections to the inverse tree-level fermion propagator 
are of C(l), while the corrections to the tree- level n — ip — -tjj vertex are suppressed by 1/A''. 



2. A further approximation concerns the self-consistent pion propagator (29). In this work four approximations for 
Gtt are considered; two local approximations and two nonlocal approximations obtained using an expansion in 
In the local approximation one parametrizes the pion propagator as 

and uses this form in all equations. In the first variant is determined as a pole mass from iG~]{pQ = M^, p = 0) = 
by the self-consistent gap equation arising from (29) 

M^=m^ + ^ {v^ + Tf{M)) - ^^fpim,) + ^^M^IpiM, 0;m,). (32) 
D •\/ N y N 

In a second variant is determined from = —iG~\{p = 0), when the gap-equation becomes 

M^ = m' + ^ {v' + Tf{M)) - ^-^fp{m,). (33) 

The subscript F denotes the finite part of the integrals defined in Eqs. (26) and (27), which are given explicitly in 
Eq. (A. 13) and Eqs. (A.16)-(A.18). In this way the finite parts of all vacuum pieces are contained in our equation. The 
importance of these terms for the thermodynamics of the PQM model was pointed out in [29] . In view of the EoS (28) 
the two definitions of coincide in the chiral limit ft. = 0, where for both variants one has = 0. We note that due 
to their self-consistent nature, when (32) or (33) is solved, a series containing all orders of I/VN is in fact resummed. 



The third, nonlocal variant of the pion equation is derived using an 1/viV expansion in the pion propagator (29) 
after exploiting the EoS (28). One obtains 

. . . =r^ + ^^^te?i + ofl). (34) 



With this form of the pion propagator the EoS reads 



+ ^ {v' + Tp{M)) + ^MM, m,) - ^^fp{m,) = ^, (35) 

where in this case = h/v and we have introduced the integral 

J(M,mg) = -i [ Gl^i{p)p^lF{p;mg). (36) 
Jp 

Solving this equation for v shows that this approximation still resums infinitely many orders in 1/ v^. 



A fourth variant of G,r j which by a strict expansion in 1 / ^/N will include terms of no higher order than 0{1/ y/N ) , can 
be obtained by expanding not only the nonlocal, momentum-dependent part of the self-energy in the pion propagator 
(29), but also its local part. This is explicitly constructed including counterterms in Sec. Ill, where a diagrammatic 
illustration of the approximation is also given. For this approximation, the pion propagator is given by Eqs. (42), (44), 
and (46), while the EoS is given in (53). 
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FIG. 1. The A dependence of the real and imaginary parts of the complex sigma pole po = Ma — iVa/2 and of the Landau ghost 
Ml in the chiral limit and for the physical pion mass. The upper curve for M^ and the lower one for represent the case of the 
physical pion mass, as indicated on the plot. Ml is shown only in this case, for in the chiral limit there is very little difference. 



E. Parametrization 



The mass parameter m^, the couplings g,X, the rcnormalization scale Mqb, the vacuum expectation value wo, and 
the external field ft-, which vanishes in the chiral limit, are determined at T = /i^ = using some information from the 
sigma sector and the following physical quantities: the pion decay constant — 93 MeV and its mass m,r = 140 MeV, 
and the constituent quark mass taken to be Mq = tun/S = 313 MeV. From the sigma sector we use the mass and 
the width of the sigma particle and the behavior of the spectral function, vq is determined from the matrix element 
of the axial vector current between the vacuum state and a one-pion state, which d\ic to the rcscaling of the vacuum 
expectation value by gives vq = /7r/2. The value of the Yukawa coupling g = 6.7 is obtained by equating the 
tree-level fermion mass rUq with Mq. The parameters A and Mqb are determined from the sigma propagator, as will 
be detailed below. Having determined them, in the chiral limit is fixed from the EoS, while in the case of the 
physical pion mass the remaining parameters and h are determined as follows. If the local approximation for the 
pion propagator is used is determined from the gap equation by requiring = m^, and h is obtained from the 
EoS. When is approximated using a large- A^/ expansion h is fixed by requiring h = m^vo, and is determined 
from the EoS. 

Now we turn to the issue of fixing A and Mqb- Using in (23) the tree-level fermion propagator together with the 
local approximation (31) for the pion propagator and also the equation of state (28), one obtains after working out the 
Dirac structure the following form for the sigma propagator: 

The integral If{p;M), obtained using the local approximation (31) for the pion propagator with = mj, can be 
found in Eqs. (10) and (11) of Ref. [59] with Mq replaced by Mqb, while ipipimq) is given in Eqs. (A.16)-(A.18). 

Both in the chiral limit M = and for M = the self-energy has along the positive real axis of the complex po 
plane two cuts above the thresholds of the pion and fermion bubble integrals, which start at p"^ = 4M^ and p^ = Arn^, 
respectively. Above these thresholds the respective pion and fermion bubble integrals have nonvanishing imaginary 
parts. Wc search for poles of the sigma propagator analytically continued between the two cuts to the second Ricmann 
sheet in the form iG~^{po = Ke~*'^,p = 0) = 0. The pole is parametrized as po = — with the real and 

imaginary parts corresponding to the mass and the half-width of the sigma particle. The solution for M„ and is 
shown in Fig. 1 both in the chiral limit [h = 111^^ = 0) and for the ft 7^ case. Similar to the case of the 0{N) model 
studied in Ref. [59] , in the chiral limit the value of is a little smaller and the value of larger than in the ft ^ 
case. Comparing Fig. 1 with Fig. 2 of Ref. [59] obtained in the 0{N) model, that is without fermions, the M„{X) 
curve moved slightly upward, but the ro-(A) curve moved significantly downward, which means that in the present case 
the phenomenologically expected value [66] Mcr/r„ ~ 1 cannot be achieved for any value of the coupling A. Another 
difference is that for low values of A there are two poles of G„ on the negative imaginary axis in contrast to only one 
such pole in the 0{N) model. These poles approach each other as A increases and after they collide at a given value 
of A there are two complex poles at higher A, one with positive and one with negative real part. The imaginary part 
of the complex pole having positive real part is shown in Fig. 1 for the rcnormalization scale Mqb = 885 MeV. As 
explained in the study done in the chiral limit in [40] for lower values of the renormalization scale the scale Ml of the 
lower Landau ghost on the imaginary axis comes even closer to M„ and as a result the spectral function of the sigma is 
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heavily distorted. In order to avoid this and based on the ratio of M„/T„ we have chosen A = 400 and Mqb = 885 MeV. 
For these values M„ = 456 MeV and = 221 MeV in the chiral case, while = 474 MeV and r„ = 152 McV for 
the case of a physical pion mass. These values are used throughout Sec. IV also in the case when the pion propagator 
is expanded in 1/y/N, that is in the third and fourth variants of discussed in Sec. II D. 



III. RENORMALIZATION 



The lesson one can learn from the successful renormalization in 2PI [67-71] or large- [64, 72, 73] approximations is 
that, due to the systematic natiu'c of these expansions, the countcrtcrms can be obtained by analyzing the structure 
of the equations. In these cases there is no need for an order-by-order detailed study of the counterterm diagrams 
which becomes rather complicated due to the proliferation of the diagrams. If one uses some ad-hoc approximation 
which spoils the self-consistent nature of the propagator equations or is not systematic, then one loses the possibility 
to explicitly or uniquely determine the counterterms from the equations. 

In this section we discuss the renormalization of the model in the case of a strict 1/VN expansion in the pion 
propagator and the equation of state. This expansion is not entirely consistent because as mentioned in Sec. II D we 
use for simplicity tree- level fcrmion propagators despite the fact that this expansion produces 0(1) corrections in the 
fermion propagator (24) which all have to be rcisummed. We will see that as a consequence of this approximation the 
best one can achieve is to determine the counterterms to some order in the Yukawa coupling. It will turn out that the 
order depends on the equation and the type of subdivergence we are looking at. 

The counterterm functional needed to renormalize the pion propagator and the equation of state reads 



N- 



5m\ 



+{N-1) — 
^ ' 24 



2 

[cAk) -{N-l)^[k'G^ 

Jk i Jk 



(38) 



Compared to the counterterm functional used in Eq. (48) of [64] to renormalize the stationary equations of the 0{N) 
model the only difference in (38) is the appearance of the term containing the wave- function renormalization counterterm 
dZ. This is needed to remove the momentum-dependent divergence of the fermionic contribution to the pion propagator 
(29) rewritten as 



iG-\k) = {l + 6Z)k'^ - M'^ 



k I{k;mq). 



The local part containing the remaining counterterms reads 

M2 = + 6ml + + ^ J G^p) - ^^r(mj. 

Using in (39) the notation introduced in (A. 15) one can readily determine 5Z : 



5Z = 



7div(fc;m,) 



(39) 



(40) 



(41) 



Separating the LO and NLO contributions in the local part given in (40) by writing = M^q + M^-^Q/y/N, an 
expansion in powers of 1/ a/TV in the pion propagator (39) gives 



G.(p) = GLo(p)-i-^'°^^^ 



N 



M^^Q-2N^g'p^lF{p;mg) 



O 



N 



(42) 



where Glo(p) = — ^to)- ^hc counterterms are also written as the sum of LO and NLO contributions Smf 



Smf°^ + 6mf^''/VN, 6Xi = 6Xf'^ + 5Xf'^/y/N and used together with (42) in (40) to obtain the equations for the LO 
and NLO local parts 



,2(1) 



,^2 2 e 2(0) A-I-^aI 2 A-l-JAn , 
Mia =vr? + ^mf + g-^« + g-^?^(^Lo), 



LO — '"- "'"•2 

1 /(0;Mlo) 



^NLO 



A 



(0) 



6 



A 



(0) 



(43a) 

i|^7VeT(mg) + ^7VeJ(MLo,m,),(43b) 



where we divided the second equation by A^^ = A -|- i5Ao°^ and used the integral introduced in (36) with replaced 
by Glo- 

Compared to the perturbative renormalization of the fermionic trace-log contribution to the effective potential per- 
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formed in [31], the difficulty here is that, due to the self-consistent nature of the pion propagator, an infinite series 
of coupling counterterms of the mesonic part, that is i5Ao, 8X2, SryiQ, and dmi has to be determined to 0{-/N) of 
the large- A'' expansion. In order to achieve this we apply the method developed in Refs. [60, 71, 73]. The method 
for determining the counterterms appearing in a particular ciquation resides in the separation of the divergent part 
of the integrals contained by the equation. This is obtained by expanding the propagators around an appropriately 
defined auxiliary propagator (sec Appendix A of [71]). Then, the finite part of the integrals is used to write down 
a finite equation, which, when subtracted from the original equation, provides a relation between counterterms and 
divergences. This relation still involves the vacuum expectation value v and the finite part of some integrals, e.g. Tp, 
the finite part of the pion tadpole. Requiring as in [73] the vanishing of the coefficient of and Tp (cancellation of 
subdivcrgcnccs) one obtains the coupling counterterms, while the remaining part of the relation mentioned above gives 
the mass counterterm (cancellation of the overall divergence). 

To apply the above method to Eq. (43a) for M^q, which is the gap equation of the 0{N) model at leading order of 
the large- TV approximation, one uses the expression of the pion tadpole given in (A. 6) in terms of finite and divergent 
pieces. Retaining in (43a) the finite part of the tadpole one obtains the LO finite gap equation 



(44) 



Subtracting this from (43a) one requires the vanishing of the coefficient of v"^ and Tf(Mlo) in the resulting equation. 
This determines the LO coupling counterterms 



A2 



T, 



(0) 



1 + Arj°V6 



(45) 



while requiring the cancellation of the remaining overall divergence determines the LO mass counterterm Sm2^^^ = 



-(A + <0 + [M' - M2]rrj /2. 

The determination of the counterterms in the equation for the NLO local part in the pion propagator parallels to 
some extent the analysis of the NLO divergences in the 0{N) model discussed in Sec. VI B of [64]. One observes that 

since /div(0;MLo) = T^^\ in view of (45) the left-hand side of (43b) is finite and it actually enters the finite equation 
for M^Lo 



A^NLO 



A 



N,JF{Mj^o,mg). 



(46) 



Subtracting (46) from (43b) the following relation between divergences and counterterms is obtained: 



A' 



1 

Toy 



dm 



2(1) , 



(1) 



+ 



,5 A, 



(1) 



g -v'' + -^T(Mlo) - 452fdiv(mg) 



4g^N,^TF{mg) + -|-^cJdiv(MLo,mg). (47) 

O D 



Then we use the expression of Jdiv(AfLO) Wg) given in (A. 28) in terms of Mloj and T(mq) in which one substitutes 
for Mlo its expression from (44). The terms proportional to Tpimq) cancel. The overall divergences determine the 
form of 5m2^^K The remaining terms proportional to v'^ and TpiMj^o) can be grouped as 



+ 



Tf(Mlo) 



+ ^ri°) + ^iV.AA-^T, 



4g*N, 

(0)^(1 



(o)^2^ 



0. 



(48) 



Requiring the vanishing of the coefficient of v'^ and Tf(Mlo) determines the NLO coupling counterterms 6X2^^ and 
^Aq^^ One can see that these counterterms agree at 0{g'^) but differ at 0{g'^). 



A completely similar analysis performed on the EoS 



m + drriQ + 



6 6 

gives an equation analogous in form and meaning to (48) 



J h 



(49) 



AW4^ (0) (0) /„(/), .„(0)x2\ R^2^(0) 
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6 



X6X 



(1) 

2 ^(0) 



21 
3 



= 0. 



(50) 



From this equation one can see that the NLO couphng counterterms 6X^2 ^ and SX^p agree at 0{g^) but differ at 
0{g*). Moreover, one sees that SX'^^ determined from the requirement to cancel the coefhcient of v'^ in (50) differs at 
0(17*) from 5X^^ needed to cancel the coefficient of Tp{Mi^q) in (50). These requirements give the same expression for 
the counterterm only at 0{g^). 

The above feature is a consequence of the fact that by keeping the fermions unresummed one does not take into 
account all the diagrams which are of the same order in the large- A'^y^ expansion. It does not necessarily mean that the 
approximation we use is unrenormalizable. Rather we suggest the interpretation that the approximation is such that 
different subseries of the counterterms are needed to cancel the subdivergences of different equations. Although we 
have a method to determine the counterterms in each equation, corrections are to be expected starting at 0{g'^). If one 
traces back the origin of the term proportional to rF(-^^Lo) in (48) one finds that it comes from the expression (44) for 
M^o used in the divergent contribution Jdiv(-WL07 ™g) as given by (A. 28). In turn, the integral J{Mi,o,mq) defined 
in (36) is generated through the expansion (42) by the one-loop fermion bubble contribution to the pion self-energy. 
But, when the first correction to the fermion propagator is included, then we have to take into account in the square 
bracket of the expanded pion propagator (42) the contribution of the two-loop self-energy 




S2(p;MLo,m,). 



This will generate in the equation for M^^q an integral similar to J{Mi^o,mq) 



(51) 



/X(MLo,m,) = -i/ / g2o(p)E2(p;Mlo,™,) 



(52) 



which is expected to have a divergence proportional to Tp{M\^o). This divergence would result in 0{g'^) corrections in 
the ^Aq"'^^ counterterm as determined from (48). Therefore, it is expected that with a resummed fermion propagator, 
as required by the large-A^/ resummation, the determined counterterms will eventually agree in all equations. It is 
rather nontrivial to check this conjecture, even at the two-loop level, because the reduction of the two-loop integral 
in (51) performed with the method of [74] produces more than a dozen scalar integrals and their contribution should 
be analyzed in (52). We have only checked that at the two- loop level the Goldstone theorem is indeed violated as 
mentioned in Sec. II D, based on the violation of the Ward identity (30). 




FIG. 2. Leading order and next-to-leading order diagrams resummed in the equation of state obtained by expanding the self- 
consistent pion propagator to first order in l/-\/iV. The tree-level pion and fermion propagators are denoted by thin and double 
lines, while the thick line represents the resummed LO pion propagator. 



Before closing this section we give in Fig. 2 the diagrammatic illustration of the equation of state obtained by a strict 
expansion to first order in 1/y/N of the self-consistent pion propagator. This corresponds to the fourth approximation 
to the pion propagator discussed in Sec. II D. Because we do not draw the counterterm diagrams, we actually obtain 
the unrenormalized EoS which reads 

(53a) 
(53b) 
(53c) 

where J(0;Mlo) = dT{M)/{dM'^), with T(M) defined in (A.6). The first set of diagrams are the 0{VN) superdaisy 
diagrams made of pions with tree-level propagators. Their resummation is clearly provided by (53b), as one can check 
iteratively. The second set of diagrams is 0(1) and contains only a single insertion of a fermion bubble. Using Feynman 



LO 



N 



V 

,,2 



A 
6 

Ag^N, 



T(Mlo)], 



NLO 



l-|/(0;AfLo) 



-T{mq) -f —J{Mi,o,mq) 
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rules, one can readily check that when the chain of pion bubbles is resummed one obtains 



/(0;MLo)T(mg) + ij(MLo,m,) 



Adding this to the fermion tadpole and dividing by v one obtains the expression of M^^o given in (53c). The finite 
version of the EoS (53) is obtained by replacing the integrals by their finite parts taken from the Appendix. 

If one would try to use the method described in this section to renormalize the EoS (28) using the local pion 
propagator (31) with a mass determined from the gap equation (32), one would encounter a subdivergence proportional 
to If{M, 0; rriq) which is not canceled. This is an artifact of the local approximation used and a self-consistent treatment 
of the propagator would unfold this into renormalizable pieces as happened already when the propagator was expanded 
consistently in 1/y/N. As we have seen, in this case only subdivergences proportional with v"^ and Tpimq) appeared. 
As a consequence, we will not attempt to explicitly construct the counterterms when using other approximate forms 
of the pion propagator given in Sec. 11 D. Fortunately, in these cases, since the propagators are of a tree-level form, 
the finite part of the integrals can be easily determined and we assume that in a given equation the subtraction of the 
infinite part of an integral can be achieved by a corresponding subset of the full series of counterterm diagrams. 



IV. THE iJ.q-T PHASE DIAGRAM 

The thermodynamics is determined by solving the field equations, i.e. the EoS (28) and the equations giving the 
dependence on T and fig of the two real mean fields $ and When the full fermion propagator is replaced by the 
tree-level one, one has in view of Eqs. (A.1)-(A.5) and (A.35)-(A.37) 



«iil-2A^.^/iv//^ 




+g^^N. 2 (r°(m,) - T,(M)) + - + 



= 0, (55) 



where = (k^ -|- m^) 2 and M satisfies either one of the gap equations (32), (33), or (43a), or the relation = h/v. 
The other equation is similar to (55), the only difference is that the derivative is taken with respect to The integral 
in (55) is the contribution of the fermionic trace-log integral defined in Eq. (A.l), while the term proportional with 
is the contribution of the quark- pion two- loop integral in (20) given in Eq. (A. 35). When solving the field equations 
for $ and we disregard for simplicity the contribution of the setting-sun and keep only the one-loop contribution 
coming from the fermionic trace-log. The complete equation (55) is solved only in one case (see the last row of Table 11) 
in order to estimate the error made by neglecting this term in all the other cases. To solve the field equations we use 
for the pion propagator a given approximation described in Sec. 11 D as will be specified below. 

The tricritical point (TCP) and the critical end point (CEP) are identified as the points along the chiral phase 
transition line of the Hq — T phase diagram where a first order phase transition turns with decreasing /Zg into a second 
order or crossover transition, respectively. In case of a crossover, the temperature of the chiral transition is defined 
as the value where the derivative dv/dT has a minimum (inflection point of v{T)), while the temperature of the 
deconflnement transition is obtained as the location of the maximum in d^/dT. The transition point in the case of a 
first order phase transition is estimated by the infiection point located between the turning points of the multivalued 
curve v{^q) obtained for a given constant temperature. Although the precise definition of the 1st order transition point 
is given by that value of the intensive parameter for which the two minima of the effective potential are degenerate, we 
adopt the definition based on the inflection point, which is also commonly used in the literature, because we are not 
computing the effective potential, but only its derivatives with respect to the fields and propagators. 



A. Phase transition in the chiral limit 



In the chiral limit we solve the EoS (28) using only the local approximation to the pion propagator (31) with 
M"^ = 0. The critical temperature of the chiral transition and the pseudocritical temperature Td of the deconfinement 
transition at vanishing chemical potential, and the location of the TCP are summarized in Table 1 for various forms of 
the Polyakov-loop potential. On one hand, one can see that with the inclusion of the Polyakov loop T-^{fJ.q = 0) and 
Ttcp increase significantly compared with the values obtained earlier in [40] without the Polyakov loop. On the other 
hand, in all cases, the inclusion of the Polyakov loop has little effect on the value of /x™^. The increase in T^iHq = 0) 
obtained with the inclusion of the Polyakov-loop effective potential is basically determined by the value of its parameter 
To, while the value of Ttcp shows no significant variation among different cases having the same value of Tq. One can 
also see, that as explained in [56], the use of the polynomial and logarithmic effective potentials for the Polyakov loop, 
that is (5) and (7), drags the value of T-^ifiq = 0) closer to the value of the parameter To than the use of ?7Fuku(^, ^) 
given in (9). In this latter case one obtains the smallest value for Ttcp- 
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193.0 
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TABLE I. The critical temperature of the chiral transition and the pseudocritical temperature Td of the deconfinement 
transition at Hq = 0, and the location of the TCP in units of MeV obtained in the chiral limit without the Polyakov loop (first 
row) and with the inclusion of the Polyakov loop using various effective potentials summarized in Sec. II B. 
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FIG. 3. Left panel: Phase diagrams obtained in the chiral limit without and with the inclusion of the Polyakov loop. The latter 
has higher Ttcp and was obtained using L'poiy(3>, with To = 270 MeV. Shown are the location of the inflection points of $(T) 
and #(T). Right panel: Chiral and deconfinement phase transitions obtained for ?7iog($, 5) with To = 208 MeV (upper curves) 
and with To(/Ltq) (lower curves). The deconfinement transition line is obtained from the infiection point of 3>(T). 



For To = 270 MeV the deconfinement transition line in the Hq — T phase diagram is above the chiral transition line in 
all three variants of the effective potential for the Polyakov loop. This is illustrated in the left panel of Fig. 3 in case of 
the polynomial effective potential, where the chiral phase diagram is compared to the one obtained without including 
the Polyakov loop. When the logarithmic effective potential J7iog($, is used either with a constant Tq = 208 MeV or 
with the /iq-dependent To proposed in [27] one finds T^ < T^ at = 0, but at a given value of the chemical potential 
the deconfinement transition line crosses the chiral transition line and remains above it for higher values of /ig. This 
is shown in the right panel of Fig. 3, where the deconfinement transition line is obtained from the inflection point of 
$(T). The transition line obtained from the inflection point of i>(T) is practically indistinguishable from the line shown 
in the figm-e. One can see that in contrast to the case of constant To, where basically the deconfinement transition line 
is not affected by the increase of /i^, with a /ig-dependent Tq the deconfinement transition line strongly bends, staying 
close to the chiral line. The two lines cross just above the TCP. 

The lowering of the deconfinement transition in the case when To(/Xg) is used and as a result the shrinking of the 
so-called quarkyonic phase was already observed in Ref. [17]. As distinguished from the mesonie phase which is confined 
and has zero quark number density and the deconfined phase which has finite quark number density, the quarkyonic 
phase is a confining state made of quarks and is characterized by a high quark number density and baryonic (three- 
quark state) thermal excitations. Based on the fact that in the PNJL model the quantity measuring the quark content 
inside thermally excited particles carrying baryon number shows a pronounced change along the chiral phase transition 
line, the region of the fj,q — T plane for which < T < Td was identified in [56] with the quarkyonic phase. The 
first numerical evidence from lattice QCD for the existence of a phase which is neither the hadronic nor the deconfined 
phase and is characterized by a high value of the quark number density was given in [75]. This could be a candidate for 
the quarkyonic phase. Further evidence for such a state was reported also in [76] within the strong-coupling expansion 
of the lattice QCD. 

Comparing our results on the phase diagram to those obtained in the chiral limit of the PNJL model one can notice 
differences of both qualitative and quantitative nature. In the nonlocal PNJL model of Ref. [24] the deconfinement 
phase transition line starts at /Ug = below the chiral transition line both for a polynomial and a logarithmic Polyakov- 
loop effective potential with To = 270 MeV, so that the two transition lines cross at finite /Xg. In our case this happens 
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only for the logarithmic potential with Tq = 208 MeV, as can be seen in Fig. 3. In [19, 24] the values of ^^(/Ug = 0) 
and Ttcp are much larger than in our case, while the value of i-i^'"^ is similar to ours. 

B. Phase transition in case of the physical pion mass 

In the case of the physical pion mass we solve the EoS (28) using each one of the four approximations to the pion 
propagator introduced in Sec. II D. Within the approximation corresponding to a strict expansion in l/\fN of the 
pion propagator and of the EoS discussed in detail in Sec. Ill there is no CEP in the Hq — T phase diagram within a 
range < /ig < 500 McV. Without inclusion of the Polyakov loop the transition at /i^ = is a very weak crossover 
characterized by a large value of the width at half maximum of —dv/dT, ~ 100 MeV. Including the Polyakov 
loop, although the width at half maximum of —dv/dT decreases by a factor of 2 compared to the case without the 
Polyakov loop, the transition remains a crossover for fiq < 500 MeV. This means that as a result of resumming in the 
pion propagator 0{1/^/N) fermionic fluctuations, while keeping the fermion propagator unresummed, the crossover 
transition at /Ug = softens and increasing jiq cannot turn the phase transition into a first order one. For the other 
three approximations, which all rcsum infinitely many orders in 1/ v^, the phase transition tiirns with increasing /ig 
from a crossover type into a first order transition and in consequence there is a CEP in the Hq — T phase diagram. For 
these cases the result is summarized in Table II for various forms of the Polyakov-loop potential reviewed in Sec. II B. 
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local, pole 


207.2 


207.7 


12.3 


39.3 


327.0 


log 


208 


local, pole 


168.0 


167.0 


*30.3 


37.9 


326.9 


log 


270 


local, p = 


209.8 


209.3 


12.1 


33.9 


329.1 


log 


208 


local, p = 


168.5 


167.0 


*42.8 


32.7 


329.0 


log 




local, p = 


168.9 


167.4 


*42.5 


25.7 


328.7 


log 


270 


large- 


209.7 


209.3 


12.0 


34.5 


329.0 


log 


208 


large- A/' 


168.5 


167.1 


*43.0 


33.0 


328.9 


Fuku 




local, pole 


191.0 


188.7 


19.8 


36.2 


326.8 


Fuku 




local, p = 


195.3 


191.2 


21.2 


31.2 


328.9 


Fuku 




large- A^ 


195.2 


191.3 


21.2 


31.8 


328.8 


poly 


208 


large- A^, full 


188.1 


183.1 


21.4 


32.2 


329.0 



TABLE II. The pseudocritical temperatures and Td of the chiral and deconfinement transitions, the half-width at half 
maximum of —dv/dT a,t jj,q — (in the cases marked with *, due to an asymmetric shape of —dv/dT the full width is 
given) and the location of the CEP in units of MeV obtained in various approximations for the pion propagator without and 
with the inclusion of the Polyakov loop. The two local approximations are defined by (32) and (33), respectively. The \axge-N 
approximation is defined by (34) and (35). Only for the result in the last row the contribution of the setting-sun was kept in 
(55). 

In the cases studied in Table II increasing /ig drives at T = the restoration of chiral symmetry via a first order 
transition at some value /x^ > Mq. Increasing the temperature /i^ decreases and the first order chiral restoration becomes 
a crossover at a much lower temperature Tcep than in the chiral case. The inclusion of the Polyakov loop increases 
significantly the value of Tcep, but as in the chiral case it has little effect on the value of /u^^^. One can see that 
neither the choice of the effective potential for the Polyakov loop nor the value of To has a significant effect on the 
value of A*g^^. Some variation can be observed among the values of Hg^^ obtained using different approximations for 
the pion propagator. The result in the last row was obtained by keeping in the field equation of the Polyakov loop 
(55) and its conjugate the contribution of the setting-sun diagram, while in all other cases only the contribution of 
the fermionic trace- log was kept. Comparing the result in the last row with that of the last row obtained using the 
polynomial Polyakov-loop potential, one can see that the error we make by neglecting the setting-sun contribution in 
all other cases is fairly small. 
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FIG. 4. Left panel: Phase diagrams obtained for the physical value of the pion mass using the local approximation to G,r with 
a mass determined by (33) without and with the inclusion of the Polyakov loop. The latter has higher Tcep and was obtained 
using ?7poiy with To = 270 MeV. The part of the phase diagram where the transition is of first order is enlarged in the inset. 

Shown are the global maxima of d<t{T)/dT and d^{T)/dT. Right panel: Chiral and deconfinement phase transitions obtained 
for [/iog($, $) with To — 208 MeV (upper curves) and with To{tJ,q) (lower curves). The deconfinement transition line is obtained 
from the global maximum of d^{T) / dT. 



The values of and Td aX — Q arc mostly influenced by the choice of the Polyakov effective potential and the 
value of To : they decrease with the decrease of Tq and by using the logarithmic potential instead of the polynomial 
one. Using the polynomial potential with Tq = 270 MeV the confinement transition line in the Hq — T plane is above 
the chiral transition line. This can be seen in the left panel of Fig. 4 where the phase diagram is compared with the 
one obtained without the inclusion of the Polyakov loop. As in the chiral case, when a logarithmic potential is used 
with either a fixed value Tq = 208 MeV or with a /Ltq-dependent Tq, the deconfinement transition line starts at /Xq = 
below the chiral one and the two lines cross at some higher value of Hq. This can be seen in the right panel of Fig. 4. 
When To(/ig) is used the two lines go together imtil they cross each other just above the location of the CEP. This 
/Xq-dependent Tq gives the lowest value of Tcep, similar to the results reported in [21] and [32]. Because of the much 
lower value of the Tcep the shrinking of the quarkyonic phase is more pronounced than in the chiral case, as the 
deconfinement transition lines approaches the fiq axis. This is even more the case here, with a physical pion mass, since 
the deconfinement transition is a crossover and as such it happens in a relatively large temperature interval. However, 
the quarkyonic phase does not vanish completely as happens in [32], where quantum fluctuations are included using 
functional renormalization group methods. 

By studying the derivatives of the v{T) and (f>(T) curves one observes in panel (a) of Fig. 5 that at low fiq it is the 
Polyakov loop which plays the driving role in the transition: for = the dv /dT is much wider and has a small peak 
in the temperature range where $(T) shows a pronounced variation. This happens only for very low values of /Xg, as 
in the region of Hq where the deconfinement transition line is a little bit further below the chiral transition line than 
for fXq = 0, such a driving role cannot be identified. For values of /Ug where the two transition line cross and also in 




FIG. 5. Evolution of the maxima of the —dv/dT and d<b/dT with the increase of the chemical potential jjbq in case of using 
J7iog($,^) with a To = 208 MeV (left figure) and with a /Uq-dependent To parameter (right figure). 
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the region where the chiral transition Hne is below the deconfinement one can see the influence of the chiral transition 
on the shape of d^/dT. This is the most pronounced in the case of the /ig-dependent Tq where the two transition hnes 
cross near the CEP. In this case the chiral phase transition plays the driving role as one can clearly see on panel (c) 
of Fig. 5 (right). In panel (d) of Fig. 5 (left) one sees that d^/dT has two peaks. In such cases, as in Ref. [28], the 
position of the higher peak is followed to determine the deconfinement transition temperature, since the first peak is a 
result of the influence of the chiral phase transition. 

From Table II one can see that there is a correlation between the strength of the chiral crossover at /ig = as 
measured by and the location of CEP: weaker crossover (larger value of F^) corresponds in general to a larger value 
of Mq^^. In the cases marked with an asterisk in Table II a dv/dT is distorted by the temperature dependence of the 
Polyakov loop, as one can see in panel (a) of Fig. 5. For this reason in this cases we denote by F^ the full width at 
half maximum of —dv/dT. In other cases one gives the half- width at half maximum. This is measured on the left of 
the maximum, because when the gap equation is used, the threshold of the fermionic bubble is generally on the right 
of the maximum and distorts the dv/dT curve. 



V. DISCUSSION AND CONCLUSIONS 

Using the tree-level fermion propagator and several approximate forms of the pion propagator obtained within a 

large- TV/ expansion, we studied in the chiral limit and for the physical value of the pion mass the influence of the 
Polyakov loop on the chiral phase transition. We obtained that only when the local part of the approximate pion 
propagator resums infinitely many orders in 1/Nf of fermionic contributions it is possible to find a CEP on the chiral 
phase transition line of the /Xg — T phase diagram. When the logarithmic form C/iog($, $) of the effective potential for 
the Polyakov loop was used with parameter Tq = 208 MeV a crossing between the chiral and deconfinement transition 
lines was observed, with the latter line starting at /ig = slightly below the former one. In this case the existence of 
the quarkyonic phase is possible. 

We have seen at the beginning of Sec. IV B that as a result of resumming in the pion propagator 0(1/ VN) fermionic 
fluctuations obtained with a strict expansion in 1/\/N, while keeping the fermion propagator unrcsummed, the phase 
transition softens. One can easily demonstrate the same feature by including the contributions of the fermion vacuum 
fluctuations and of the pion tadpole in the equation of state of Ref. [27], that is the field equation determining the 
chiral order parameter. There, bccaiisc the parameters of the PQM model were determined at trc!(;-k!vc;l. the fermionic 
vacuum fiuctuations coming from the fermion tadpole (T°) were neglected, while the pions were treated at tree-level. 
However, by choosing an appropriate renormalization scale one can arrange for the vanishing of the entire Tp only 
at T = 0. At finite temperature the vacuum fluctuation is in this way correctly included and due to the temperature 
dependent fermionic mass Tp will be nonvanishing. The value of the renormalization scales for which the fermion and 
pion tadpoles vanish at T = ;Ug = are Mqf = y/eniq and Mqb = \/em,r, respectively. The importance of including the 
vacuum fluctuations was discussed also in [29, 31], where the effect on the location of the CEP and on the isentropic 
trajectories in the Hq — T plane was shown. From Table III one can see that comparing with the original result of 
Ref. [27] the inclusion of the fermionic vacuum fluctuations softens the transition at fXq = 0, as shown by the larger full 
width F^ at half maximum of —dv/dT, and in consequence the location of the CEP is moved to higher values of fiq 
and lower values of T. Inclusion of the pion vacuum and thermal fluctuations in the equation of state through a pion 
tadpole further accentuates this behavior. Inclusion of the fluctuations using functional renormalization group methods 
also pushed the location of the CEP to higher values of fiq, as can be seen by comparing the left panel of Fig. 6 in [32] 
to Fig. 6 of [27]. 

















TCEP 




QP 






+ 




184.6 


4.6 


162.8 


165.1 


QP 






+ 


+ 


180.2 


8.6 


145.3 


204.3 


QFT 


+ 




+ 




173.0 


26.9 


91.3 


241.1 


QFT 


+ 


+ 


+ 


+ 


170.1 


30.3 


85.5 


243.5 



TABLE III. The pseudocritical temperatures of the chiral transition and the full width at half maximum of —dv/{dT) at 
fiq = 0, and the location of the CEP in units of MeV in various treatments of the model with a physical pion mass. The Polyakov 
loop is included using f7poiy($,$) and To = 208 MeV. QP stands for the quasiparticle approximation in which the vacuum 
fluctuations in the fermion [Tp) or pion (Tp) tadpoles are disregarded (marked by — ) and only the finite temperature part of 
the tadpoles (T*^ or T^) is kept (marked by -|-). QFT stands for a quantum field theoretical calculation where the vacuum 
fluctuations are properly treated. The first row is the reproduced result of [27]. 

It remains to be seen to what extent our results are stable against the use of the self-consistent propagator for 
fermions, as required by a completely systematic large-A^/ expansion. A highly interesting question which requires 
going beyond the level of approximations of this work is whether a completely systematic expansion in 1/ ^/N of the 
propagator equations could lead to the existence of the CEP in the phase diagram, and how the results obtained 
within such a resummation scheme are related to a numerically even more demanding resummation represented by 
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the complete self-consistent solution (without further expansion in of the coupled pion and fermion propagator 

equations. 
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Appendix: Integrals and performing the color trace 
1. The trace-log 

We calculate first an integral appearing in (18) when the approximation G{k) — >■ D{k) is used. This is defined as 

h = ^/Nitrn,c j\^D-\k) = -2s/NTY,Y. J ^'4^^' H + ^1)] , (A.l) 



=1 n 



where wc used (21) and the notation i?^ = k^ + rn^. The Matsubara frequencies are ojn = {In + 1)txT — ifM, with the 
color-dependent chemical potentials defined in (15). Doing the Matsubara sum and an integration by parts one obtains 



(A.2) 



where ,f^{E) = l/(exp[/3(£^ =F /Xj)] + 1). Using the diagonal form of the Polyakov-loop operator given in (2) one can 
define following Ref. [16] 



1 ^= 1 1 

1 1 1 

i=l 



Then, simple algebra shows that upon working out the traces f^{E) can be expressed in terms of $ = (trcL)/A^c and 
l> = {ti cL'^)/Nc as 

f+(E) = !— — CA 4a) 

' l + 3(i + $e-'3(-E-M,))e-/3(^-''5)+e-3^(^-''o)' ^ ' 



Through these functions the integral Ii in (A.l) is also expressed in terms of $ and $ as 



7. = -2N^Vn j E, + — (UiE,) + /, (E.)) 



(A.4b) 



(A.5) 



(27r)3 

The derivative of this integral with respect to $ is used in (55). 



2. Tadpole and bubble integrals 

The pion tadpole integral is given by 

T(M) = / G^,i{k) =tY,I = TMM) + Tf{M), (A.6) 

where w„ = 2nnT and El = k"^ + M^. Depending on which one of the approximated pion propagators discussed in 
Sec. II D is used, the T-dependent mass M satisfies either one of the gap equations (32), (33), or (43a), or the relation 
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= h/v. Using a 4D cutoff A, the quadratic and logarithmic divergences, 

1 , A2 



7^(2) ^ _}_ 
<^ 167r2 



are separated in the divergent part 

Tdiv (M) = T^') + - M^]T^°'> , (A.8) 
by expanding, as in Appendix A of [71], the propagator Gj^j{k) around the auxihary propagator 

Goip) = (A.9) 

Here, Mq = M^^/e, where Mqb is the renormahzation scale introduced in Sec. HE. 

The finite part of the tadpole is written as a sum of two terms, having zero and one statistical factor, respectively 

Tf{M) = T^{M) + tI^{M), 

where f{E) = l/{exp{l3E) - 1). 



The fermion tadpole integral defined in Eq. (26) is written as 

N,. 

^K) = ^ E E / 1??T772 = ^^-("""^ + ^^■'^^^ 

where i?^ = + m^. ix)„, = (2n + l)7rT — i/ij. The sum over the color degrees of freedom was done with the help of 
(A. 3). The divergent part of T(mg) reads 

TdivK) = ri^) + [ml - Ml\Tf\ (A.12) 
while the finite part is decomposed as in (A. 10): 

fir(m,) = f'^(mg) + f^{mq), 

with fiiEk) given in (A.4). 



The fermion bubble integral given in Eq. (27) is decomposed as 

where i^f = k^ + m^, £"1 = (k — p)^ + m^, and a;„ = (2n + l)7rT — i/ij. The divergent part obtained after doing the 
Matsubara sum is independent of the momentum and is given by 

/div(p; m,) = Tf\ (A.15) 

while the finite part is written again as a sum of two terms, having zero and one statistical factors, respectively: 

If{p; mq) = i%{p; nig) + F{p] lUq). (A.16) 
The finite part with no statistical factor is given by 

9^ < 4m^, 



1 TO^ O f 2arctanQ \ Q < < Ami, 



Q + 1 



where Q = y^jl — 4m^/p2|, while, after summing over color indices, the term with the statistical factor reads 



1 



k 4:ElE2 



Po - El - E2 + ie po + Ei+ E2 + ie po - Ei + E2 



Po + El — E2 
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.(A.18) 



3. The setting-sun integral 



Using the definitions of the fermion bubble integral given in (27) together with the relations (A. 14) and (A. 15), the 
integral defined in (36) reads in terms of the propagator Dq introduced in (25) as 



^ (s{M,mq)-T{M)Tl 



1 + 



(A.19) 



Here, we have used that —i ;(p) = dT [M) / {dM"^) , with M satisfying the relation = h/v or the gap-equation 

(43a) and we introduced the setting-sun integral 



Nc r 

S{M,mq) = ^Y. [ [ Doik)G^Aq)Do{k + q) 



(A.20) 



With the method described in Refs. [70, 78] the setting-sun integral can be decomposed as a sum of terms containing 
zero, one, and two statistical factors: 



5(M, ruq) = S°{M, rUg) + S^'\M, rrig) + S^'^{M, rUg), 



with 



S°{M,mg) 



(0) .(0) 



Do{k)G^Aq)Do{k + q), 



(A.21) 



(A.22a) 



N, 



k J q 

J2 P{k;mq)aik)-2I^^Ak)ai{k) = T^(M) + 2r^(m,) T^^^ + ^^'^(M, m^), (A.22b) 



E 

i=l 



c „(0) .(0) 



k J q 



[d,{k)a,iq)G„,l{k + q)- 2di{k)a{q)Do{k + q)] , 



(A.22c) 



where di{k) = e{ko)po{k)[f^ {\k()\) + /- (|fco|))/2 and a{k) = e{ko)po{k) f {\ko\) , with e(fco) the sign function, and 
Po{k) = 27re(fco)(5(fc^ — m^) and po{k) — 27re(fco)(5(fc^ — M^) the free spectral functions. Here, we have introduced the 

notation j^^^ = J and we have separated in 5^'^ the divergence coming from (A. 14) and 



(0) 



lH^ik) = -I I G.,iiq)Doik + q) = Tf + /0,,^(fc), 



(A.23) 



where the finite part is 



IqiT,F{k) 



327r2 



In 



fc2 ) 



+ 



G 



167r2/!c2 



1 Af ^ + w?' -\- G 

2 1" M2 + „,2 _ fc2 „ g - - + fc2 > (M + mqf, P<{M- mq)\ 



k- + A//(2 p - Am2 

arctan — \- arctan - 



(A.24) 



(M - mq)2 < fc2 < (M + mg)2. 



G G ' 

with Am2 = M2 - m2 and G = y/p^ _ 2^2(^2 + ml) + (M2 - m2)2. 

In S"{AI,m.q) one expands, as in Ref. [71], the propagators G^r,; and Dq around the auxiliary propagator defined in 
(A. 9) and one obtains 



S°{M, ruq) = S°{Mo) + T\M) + 2f\mq) - ST^ tJ°^ + [M^ - + 2{ml - M^)] + 5° (M, ruq), (A.25) 



2M rr{i) I cO , 
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where 



fiO) 



AO) 

I Go{q)Goiq+p)-T^ 

Jq 

The divergent parts of 5° and S^'^ combine, so that the complete divergence of the setting-sun integral reads 

5div(M,mg) = 5°i,(M,mq) + S^(^{M,m,) = S°{Mo) + [t{M) + 2f{m,) - 3Tf ] Tj°) + [M^ + 2ml ' 3^o)] T^'^ ■ 
Using this expression in (A. 19) one obtains 

Jdiv(M, m,) = S°{Mo) + \2f{mg) - 3Tf ^1 T^°^ + [2{M^ + m^) - 3M^]T^'^ 



3 dM§ 



167r2 



(A.26) 



(A.27) 
(A.28) 



In what follows wc give the finite part of the setting-sun integral in terms of which Jf{M, rUq) can be easily obtained. 
With the help of (A.26) one obtains from (A. 25) 



S%{M, mg) = 5°(M, rUg) - 5°(Mo) - ^ [M^ + 2ml ' ^^o] 



(0) 



167r2 



ivr In 



(A.29) 

Sp{M,mq) can be most easily calculated numerically as in [77] by going to Euclidean space and using 4D cutoff A to 
regularize the integrals. By doing the angular integration one obtains 



1 r^" r^" ' 

S\M,m,) = -^^ dxj^ dy- 



■y + m^ - ^{x + y + mlY - 4xy 
{x + ml){y + M^) ' 



(A.30) 



with a similar integral for S^{M), but with replaced by Mq. Then, by calculating tJ°^ defined in (A. 7) with the 
same cutoff A, one can look for the range of A where S%{M,mq) is insensitive to the variation of A. In this work we 
have used A G [600, 800] GeV. 

The finite part of the setting-sun integral with one statistical factor factorizes upon integration over fco due to the 
Dirac delta's of the free spectral functions and one has 



S^/{M,mg) = J° (fc2 = M^,mg)T^{M) + 2I°^^p{e = ml)f^{mg). 



(A.31) 



After performing the frequency and angular integrals in the part of the setting-sun integral containing two statistical 
factors one obtains 



2 



1 

E1E2 



[f+{E,) + fr{E,))f{E) 



In 



{El + Ef - £2 



{El + Ef - E^_ 



In 



{El - Ef - El 



{El - Ef - Ei 



{j+{Ei)fr{E2) + fr{Ei)f+{E,)) In 
'f+{Ei)f+{E2) + fr{Ei)fr{E2)) In 



{Ei+E2f-^ 



{El + E2f - E. 

{El - E2f - El 



{El - E2f - Et 



(A.32) 



where Ei = ^k^ + m^, E2 = ^cf + m^, E = Vq' + M\ and E± = ^/{\k\±\^fTW , E± = ^(|k| ± |q|)2 + m^. 

For the terms containing one fermionic statistical factor the sum over color indices will give the functions introduced 
in (A. 4). For the terms containing two fermionic statistical factors one introduces the notation X± = exp[(3{Ei =F /x^)] 
and Y± = exp[/3(£'2 T Mq)] and writes the statistical factors in terms of the Polyakov-loop operator. One obtains 



J2ftiEi)fr{E2) = ^J2'r^ 



1 



1 



LX+ + 1 L^Y_ + 1 



i=i " i=i 

1 + x+y_(i + X+Y_) + $ [2y_(l + X+Y_) +xl]+^ [2X+(1 + X+) + F?] + 3$$x+y_ 



(A.33) 



{Xl + 3^X1 + 3^X+ + 1){Y^ + 3$F2 + 3$y_ + 1) 
For the term with ff{Ei)f^~{E2) one obtains the same expression, but with X+ replaced by y+ and y_ replaced by 
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X_. Similarly, one has 
1 



i=l 



3$2j^2y2 ^ 6$2x+F+ + (3$$ - 1)X+Y+{X+ + y+) + 2$(X+ + Y+- XlY^) + - Y+f + 1 

(X3 + 3^X1 + 3^X+ + l){Yf+ 3$F2 + 3$y^ + 1) 



(A.34) 



and for the term with {Ei) f^ {E2) one has the same expression, but with $ interchanged with and X+ replaced 
by X- and y+ replaced by Y-. We note here that, as one can see by setting $ = $ = 0, the expression on the right- 
hand side of (A.34) can be negative and in consequence it does not allow for an interpretation in terms of distribution 
functions. 

When the approximation G{k) D{k) is used, the setting-sun integral defined in (A. 20) appears also in the expression 
of the quark-pion two-loop integral of (18) which reads 

/2 = _V]v^itrc,c / / 75-D(fc)75£'(fc +p)G,r,i(p) = sWat, ^(m,) - 2T(mg)T(M) - M^S{M,mA , (A.35) 



k J p 



where depending on the approximation on the pion propagator M satisfies either one of the gap equations (32), (33), 
(43a) or the relation = h/v. Here, we introduced the integral 



which, after performing the Matsubara sums, can be written as 



Tf^(m,) 



167r4 7o 



/•oo ,22 



(A.36) 



(A.37a) 
(A.37b) 

(A.37c) 



where Ei = yjk'^ + rri^, E2 = y^q^ -|- m^, and the color sums are given in (A. 33) and (A.34). In terms of the quantities 
defined in (A. 37) the derivative with respect to $ of I2 appears in (55). 
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